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Abstract. We study dynamics of A(t) cosmological models which are a natural generalization 

of the standard cosmological model (the ACDM model). We consider a class of models: the 

2 

ones with a prescribed form of A (t) = Abare + jt- This type of a A(f) parametrization 
is motivated by different cosmological approaches. We interpret the model with running 
Lambda (A(t)) as a special model of an interacting cosmology with the interaction term 
— dA(t)/dt in which energy transfer is between dark matter and dark energy sectors. For the 
A(t) cosmology with a prescribed form of A(f) we have found the exact solution in the form 
of Bessel functions. Our model shows that fractional density of dark energy U e is constant 
and close to zero during the early evolution of the universe. 

We have also constrained the model parameters for this class of models using the as¬ 
tronomical data such as SNIa data, BAO, CMB, measurements of H(z) and the Alcock- 
Paczynski test. In this context we formulate a simple criterion of variability of A with 
respect to t in terms of variability of the jerk or sign of estimator (1 — O mj o — The 

case study of our model enable us to find an upper limit a 2 < 0.012 (2cr C.L.) describing 
the variation from the cosmological constant while the LCDM model seems to be consistent 
with various data. 
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1 Introduction 


The standard cosmological model describes the matter content of the Universe comprising 
the cold dust matter (baryonic matter and dark matter) which satisfies the equation of 
state for dust p = 0. In turn, dark energy is described in terms of an effective parameter 
(the cosmological constant) which should be treated as the best ‘economical’ (as only one 
parameter used to describe the whole dark sector) description of the cause that the Universe 
expansion accelerates in the current epoch. 

The natural interpretation of the cosmological constant arises as an effect of quantum 
vacuum energy. Since this form of energy should be independent of the reference frame it must 
be proportional to the only ‘invariant’ second order metric tensor i.e. = p vac <7/w If 
we include the conservation condition which for the cosmological model with the Robertson- 
Walker (R-W) symmetry assumes the form 

P = ~3H(p + p), (1.1) 

then we obtain that p vac = const = A and p VSLC = —A, H = ^(ln a) is the Hubble parameter, 
where p is total energy density, p is total pressure, an overdot denotes differentiation with 
respect to the cosmological time f: we use a natural system of units in which 8irG = c = 1. 

If we interpret the cosmological constant A as a vacuum energy, then there is a difference 
between its today value required to explain observations of type la supernovae (SNIa) and 
the value of p v ac estimated from effective field theory. The former is smaller by a factor of 
1CU 120 . This discrepancy is called the cosmological constant problem. 

To achieve the conservation of energy-momentum tensor (divergence of energy-momentum 
tensor T^ v is vanishing) different descriptions of dark energy sector have been proposed. In 
the simplest case the time cosmological term A it) is shifted to the right-hand side and treated 
as a source of gravity. Such an approach is called a A(t)CDM cosmology. 

In this paper we assume p = p m + pd e where p m is a density of matter and pd e is 
the density of dark energy. We also assume p m = Pb + Pdmi and p m = Pb + Pdm-, where 
Pb = Pb,oa(t)~ 3 and pb = 0 are a density and a pressure of baryonic matter, pd m and pdm = 0 
are a density and a pressure of dark matter, a(t) is the scale factor. The state equation 
for dark energy is assumed as pd e = — Pde- I n this case the conservation condition has the 
following form 

Pdm + 3-ffpdm = Qi (1-2) 

Pb + 3Hp h = 0 (1.3) 
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and 


Pde — Qi ( 1 - 4 ) 

where Q describes an interaction between dark matter and dark energy and this case is 
expressed by Q = —A. The conservation condition can be rewritten in the form 

Pm = -3 Hp m - A, (1.5) 

where p m and p m are energy density and pressure of matter. 

Pani et al. considered the energy-momentum tensor which ensures the covariantness 
of general relativity [1]. An alternative approach is to postulate the scalar field (j> with the 
potential V(cf>) for this model [2] which guarantee that model is covariance. 

We consider a model with a parametrization of A following the rule 

a 2 

A (t) = Ab are + (1-6) 

where a 2 is a real constant; Abare is a constant and p mc = A. This model belongs to a 
larger class of cosmological models with interaction. In this case the interaction term is 
Q = —dA/dt. In this model the interaction is between dark matter and dark energy. This 
model belongs to a class of models so-called early constant dark energy during the matter 
dominating stage. 

If we replace the cosmological time t by the Hubble scale time in eq. (1.6), then we 
obtain the A (H) parametrization which is based on Lima at al. [3-5]. 

We estimate the value of the parameter a 2 as well as the other models parameters from 
available astronomical data. This class of models is compared with the standard cosmological 
model (the ACDM model). 

Let us enumerate motivations for introducing form (1.6) of parametrization of dark 
energy. 

1. The parametrization of dark energy can be derived from the quantum mechanics which 
describes how decaying false vacuum states changes in time. It can be shown that at the 
late time it can be identified as the cosmological constant which is time dependent and 
changes following the rule (1.6) and parameter a 2 is small and constitutes a leading 
term for long-term behaviour in power series of energy density of decaying vacuum 
[2, 6][7, 8], 

2. A new model of agegraphic dark energy [9, 10] based on some quantum arguments that 
the energy density of metric fluctuation of the Minkowski spacetime is proportional to 
h and and it also motivated Karolyhazy uncertainty relation [11]. If we identify the 
time scale as the age of the Universe T, then we obtain that the agegraphic dark energy 
is p q oc ^r. 

3. In the de Sitter universe there is a possibility to define in the framework of general 

relativity length and time scales A (t) = = c 2 p ^ [12]. Otherwise, any cosmo¬ 

logical length scale or time scale can determined the relation A(t). Chen et al. [12] 
demonstrated how holographic [13, 14] and agegraphic dark energy conceptions can 
be unified in the framework of interacting cosmology in which the interacting term is 
Q = — p\. The variational approach to an interacting quintessence model was recently 
considered by Bohmer et al. [15]. 
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4. Ringermacher and Mead [16] considered dark matter as a perfect fluid satisfying the 
equation of state p = — \p. The energy density of such fluid mimicking dark matter 
effects varies like A rather than d | rk as in the standard cosmological model. 

5. Haba has discussed recently cosmological models of general relativity in which a source 
of gravity (right-hand sides of the Einstein equations) is a sum of the energy-momentum 
of particles and the cosmological term describing a dissipation of energy-momentum. 
He obtained a cosmological model with the cosmological term decaying as 1/f 2 [17, 18]. 


2 Exact solutions for A(f)CDM cosmology with A (t) = A + p- 

For the parametrization of A (t) (1.6) it is possible to obtain exact solutions and discuss 
cosmological implications of this generalized standard cosmological model. We show that a 

deviation of this model from the ACDM model can be probed by a measurement of a jerk. 

2 

We start from the Friedmann first integral in the FRW cosmology with A(f) = Abare+^j 
where t is the cosmological time and a 2 is either positive or negative, 


a 


3 H(t) — p m (t ) + Abare + 


and the conservation condition 


Pm(t ) = -3 H(t)p m (t) - 


d( Abare + Jl) 

dt 


Equation (2.1) can be rewritten in the dimensionless parameters 


f^m,0 — 


Pm,0 

3 W 


^A,0 — 


Abare 

W 




a 


zh 2 t r 


( 2 . 1 ) 


( 2 . 2 ) 


(2.3) 


where To is the present age of the Universe, i.e. To = J 0 T ° dt = . and quantities labeled 

by index ‘O’ are defined at the present epoch for which ao = 1. Then equation (2.1) has the 
following form 

H{ty = Q m (t) + U A , 0 + 1*^, ( 2 - 4 ) 


H 2 


where U m (t) = flb.o a (t) 3 + fldm,o f(t) and /(To) = 1. At present, equation (2.4) is expressed 
by 

1 = flm.O + ^A,0 + ^o,0- 

After differentiation of both sides of (2.1) with respect to t we obtain 


(2.5) 


6 H(t)H(t) = p m {t) + 


d{ Abare + jr) 
dt 


Equation (2.6) can be simplified with the help of (2.2). Then we obtain 

After substitution of (2.1) to (2.7) we obtain 


a 


H{t) — — ( Abare + ^- 3H(t ) 


( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 
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Equation (2.8) can be rewritten in the dimensionless parameters. Then we obtain 


h(t) = 3 ^(n A, 0 + ^ 2 — 2 


- Kty , 


(2.9) 


where h(t) = j 

The general solution of equation (2.9) has the following form 


, , . 2 d , 

m = sw a dt hg 


'Mx .(vs; 


( 2 . 10 ) 


where C\ is a constant, an J n (x) and Y n (x) are Bessel functions of the first and second kind, 
the index n of these functions is given in terms of H a> o, and To, n = \\/l + 9Q a fiTQ . 
We can rewrite (2.10) to the form 


, , , 2 d, 

'* <t) = 3^* log 


St Ur. + n 2 j„ VEW 


( 2 . 11 ) 


For the correspondence with the ACDM model (a 2 = 0) we choose D\ = 0. Then solution 
(2.11) is given by the formula 


w , 2 d, 

^ = 3H 0 dt'° g 


StliJ hYYEi, 


( 2 . 12 ) 


where I n {x) is the modified Bessel function. Solution (2.12) can be rewritten to the following 
form 

j ( Sy/flAfiHo 
1 n— 1 1 9 ‘ ) 

(2.13) 


1 _ O n ,_ 1 \ 2 

h(t) = OZJ + + \/^A,cr 


3Hot 


In 


3y/n A , 0 H 0 


Because H(t ) = 4? In a, then it is easy to obtain the scale factor from (2.12) in the form 


a(t) = C 2 


St 


(2.14) 


The diagram of a(t) is presented in figure 1. We obtain a formula for p m (t ) from (2.7), (2.9) 
and (2.12) 


Pm — —3 Hq I 11a,0 + 


Ho,oTo 


2 d 
3 Hq dt 


log 


st(u 


(2.15) 


or in an equivalent form 

( 


Pm ~ ~3 Hq 


r> | ^, 0 T 0 2 

n Y0 + t 2 


V 


1 - 2 n 
3H()t 


+ \/ ^A,0 


In —1 




T ( 3y/£lA,oHa + 

V 2 J ) j 


(2.16) 
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The diagrams of p m (t) and p m (a) are presented in figure 2 and 3. In comparison, the diagram 
of pde(t) is demonstrated in figure 4. 

The dark matter is expressed by 

Pdm = Pm ~ Pb,oa~ 3 - (2.17) 

If we use formulas (2.14) and (2.16) in equation (2.17) then we get 


Pdm 


= -3//,) 


^A,0 + 


^ 0 , 0^0 


1 - 2 n 
3Hot 


+ 


In-, 


- PbfiC 2 


1 —3 


T ( 3y/fl A ,oHo ^ i 

n v 2 ) J ) 


1 -2 


. (2.18) 


The ACDM model can be obtained in the limit fi a ,o = 0. Then index n = \ and 
Ii (x) = J ^ skill (x). Finally the solution (2.12) reduces to 


, , . 2 d, 

h{t) = W 0 dt log 


sinh 




Equation (2.19) can be rewritten to the equivalent form 


h{t) = v^a.o coth 




In the special case the solution of (2.9) for 12 a, o = 0 has the following form 


hit) = 


3 Ht 


1+ /l + 9ij2o qT 2 1 _ 


Ci 




(2.19) 


( 2 . 20 ) 


( 2 . 21 ) 


For the correspondence with the CDM model we choose C\ = 0. Then equation (2.21) is 
simplified to 

' ( 2 . 22 ) 


h{t) = 


3Hot 


l + ^/l + 9H$n Qfi T$ 


From equation (2.22) we can obtain an expression for the scale factor 

a(t) = c 2 t*( 1+y / 1+9H ° na ’ oT °). 


(2.23) 


If we know an exact solution for the scale factor a(t) it will be possible to calculate a 
dimensionless parameter called a jerk related with a third order time derivative of the scale 
factor 


3 = H(t) 3 a(t) [~dt 3“J ' (2 ‘ 24) 

After some calculations we obtain the third order time derivative of the scale factor in the 
form 


a = 


3 H$a 


n A0 + ^201 - h A + :iH JC _ 2 hh) + H 2 ah 2 + 2H 2 ahh. (2.25) 


3Hh 


212„ nTn 
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a(t) 



Figure 1. Diagram of the scale factor a[t) for three cases. The top function is for a 2 = —0.2, the 
middle function represents the ACDM model and the bottom function is for a 2 = 0.2. We assume 
that Hq = 68.27 km/(s Mpc) and Q rn = 0.35. Time t is expressed in a unit (100 s Mpc/km). 


Pm(X) 



Figure 2. Diagram of the energy density p m (t) for three cases. The top function is for a 2 = 0.2, the 
middle function represents the ACDM model and the bottom function is for a 2 = —0.2. We assume 
that Hq = 68.27km/(s Mpc) and D mj0 = 0.35. Time t is expressed in unit (100 s Mpc/km). We 
assume 87 tG = 1 and we choose for p m a unit (km/(100 s Mpc)) 2 . 


A substitution of the expression h(t) from (2.9) gives us the exact formula for the jerk as a 
function of the cosmological time t 


j{t) = 1 - 


H 0 P 


( 

1 - 2 n 
3Hot 

V 


In- 


+ 


A,0“ 


(3y/Cl A O H 0 \ \ 

1 v— 5 —V 


-3 




(2.26) 
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Figure 3. Diagram of the energy density p m (a) for three cases. The top function is for a 2 = —0.2, 
the middle function represents the ACDM model and the bottom function is for a 2 = 0.2. We assume 
that Hq = 68.27km/(s Mpc) and £l m ,o = 0.35. We assume 87 tG = 1 and we choose for p m a unit 
(km/(100 s Mpc)) 2 . 


Pde(t) 



Figure 4. Diagram of the energy density Pde{t) for three cases. The top function is for a 2 = 0.05, 
the middle function represents the ACDM model and the bottom function is for a 2 = —0.05. We 
assume that H 0 = 68.27km/(s Mpc) and fl m ,o = 0.35. Time t is expressed in unit (100 s Mpc/km). 
We assume 87 tG = 1 and we choose for pde a unit (km/(100 s Mpc)) 2 . 


The jerk calculated for t = To, i.e. for the present epoch is given by formula 


Jo = 1 - 


HoTo' 


(2.27) 


The diagram of j(z) is presented in figure 5. One can see the jerk can be treated as a tool 
for detection the variability of dark energy. 

From the exact solution (2.26) one can see that the deviation of the generalized model 
from the ACDM model is given by time dependent contribution to the jerk because for the 
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j(z) 



Figure 5. Diagram of j(z) for Hq = 68.27 km/(s Mpc) and f2 mj o = 0.35. The top function is 
for a 2 = —0.05, the middle function represents the ACDM model and the bottom function is for 
a 2 = 0.05. 


Om(z) 



Figure 6. Diagram of Om(z) for Hq = 68.27 km/(s Mpc) and D m ,o = 0.35. The top function is 
for a 2 = —0.05, the middle function represents the ACDM model and the bottom function is for 
a 2 = 0.05. 


ACDM model it is equal one. Therefore if we can detect from the astronomical observations 
the time variability of the jerk it will be a simple diagnostic of decaying vacuum. If D Qj o is 
non-zero this means that D mj o + Da,o < 1- Note that 

a 2 

O*,0 = 1 — ^m,0 — ^A,0 = g fl2rp2 ' (2.28) 

Because To < i.e. < 1 and a 2 = 3//qTqD Q] o, i.e. 

a 2 

y < (2.29) 
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2 

and from the estimation of one can obtain an upper limit on 

From the value of the jerk for the current epoch (see formula (2.27)) there comes also 
the limit values of the jerk 


1-3D Qj0 < j 0 < 1, 

for a 2 > 0 

(2.30) 

1 < jo < 1 - 3O q , 0 , 

for a 2 < 0. 

(2.31) 


Sahni et al. [19-21] proposed in the context of testing and comparison of alternatives 
for the ACDM model Om(z) diagnostic test 


Om(z ) 


h 2 (x ) — 1 
x 3 — 1 


(2.32) 


where x = 1 + z. While this parameter is constant for the ACDM model, Om(x) = D nij o for 
any deviation from zero would discard the ACDM model for the description of the cosmic 
evolution of the current Universe for low z. But Om(*) diagnostic test is not constant for high 
z because ACDM model should respect radiation for high z. Note that if the radiation density 
is included then the behavior of Om(z) will be different for the case of matter and cosmological 
constant. For high redshift the contribution from radiation density will dominate. In our 
paper matter and energy density is present at very beginning and effect of radiation density 
is not included because of complexity of analytical calculations. Therefore our comparison 
of a jerk and Om(z) is not valid for high redshift. Let us note that in our case Om(x) is not 
constant and evolves with the cosmological time as 


Om(t) 




In— 1 




7 


-1 


vtJiJ 


’Tr 


2 pt (/„ 



(2.33) 


The diagram of Om(z) is presented in figure 6. From comparison of figures 5 and 6 one can 
observe two alternative ways of the detection of the variability of dark energy with respect 
to time. 


3 Dynamics of the generalized ACDM model 

For a deeper analysis of dynamics it is useful to investigate how exact solutions (trajectories) 
depend on initial conditions. The natural language for such a discussion is the phase space 
which a space of all solutions for all admissible initial conditions. 

Let us consider now the dynamics of the model under consideration as a dynamical 
system. In this paper we consider the case of a positive cosmological constant Abare > 0 
and strictly positive energy density of matter p m > 0. The first step in a formulation of the 
dynamics in terms of a dynamical system is a choice of the state variables. Assume state 
variables are as follows 


x 2 = 


Pm 

37P’ 


y 2 = 


A 


bare 


3H 2 


z 2 = 
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3 H 2 t 2 ' 


(3.1) 




















We also choose a new time variable r: r = In a; let a prime denotes the differentiation with 
respect to the Hubble time r. Then, we differentiate with respect to r the expressions for 
x 2 , y 2 and z 2 in (3.1) and obtain 


2xx' = 


2xx 

~JT 


Pm 


0 / 2 yy 

2 » = it = ~ 


_ 2pmH 

3 H 3 3AT 4 ’ 

2Abare^f 


3 H 4 
2 


2zz > = ^ = _ 

H 3H 3 t 3 


2 H 

3 HH 2 ' 


(3.2) 

(3.3) 

(3.4) 


Due to relation (2.2) the expression for p m can be replaced by —3 Hp m + and then with 
the help of (2.7) H can be replaced by — As a consequence we obtain the set of equations 


„ / _ Pm 2a 

2XX ~ H 2 + 3H 3 t 3 


+ 


2 Pm 
6/7 4 ’ 


2y?/ , = 


2Ab are p m 


2zz' = — 


6i7 4 

2 


+ 


2pn 


3H 3 t 3 ' 6HH 2 ' 

After returning to the original variables x, y , z we obtain the system 


(3.5) 

(3.6) 

(3.7) 


/ 3 r 2 £ 3 3 o 

X = --x + V3a-h -3T, 

2 x 2 

/ 3 2 

y = 2 X y> 

z' = —-s/3 z 2 + ^x 2 z. 


(3.8) 

(3.9) 
(3.10) 


Note that the right-hand side of (3.8) is not defined on the plane x = 0. All state variables 
are constrained by the condition x 2 + y 2 + z 2 = 1, i.e. phase space is a surface of a three- 
dimensional sphere. 

We regularize system (3.8)-(3.10) in such a way that its right-hand sides are in a polyno¬ 
mial form. For this purpose we introduce new state variables X,Y, Z: X = x 2 ,Y = y, Z = z. 
Note that transformation x —> X is not a diffeomorphism on the line x = 0. Then system 
(3.8)-(3.10) represents the dynamical system with smooth right-hand side functions, namely 


X' = -3X + 3X 2 + 2/3 a 2 Z 3 , 



Z' = -y/3 Z 2 + ^ZX, 

where the phase space is restricted by the condition 

X + Y 2 + a 2 Z 2 = 1. 


(3.11) 

(3.12) 

(3.13) 


(3.14) 


The critical points of the system (3.11)-(3.13) , their type and dominant contribution in the 
energy constraint X + Y 2 + a 2 Z 2 = 1 are presented in table 1. 
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Table 1. The critical points of the system (3.11)-(3.13), their type and dominant contribution in the 
energy constraint X + Y 2 + a 2 Z 2 = 1. 


No 

position of critical point 

type 

dominant contribution in (3.14) 

1 

Xq = 0, Yq = 1, Zo = 0 

stable node 

A dominant state in 
the future (de Sitter) 

2 

X() = 1, Yq = 0, Z 0 = 0 

unstable node 

matter dominant state in 
the past (Einstein-de Sitter) 

3 

= 3^2 ( — 1 + Vl + 3a 2 ) 

To = 0 

Z ° — V 3 a 2 ( 1 + Vl + 3a ) 

saddle 

both decaying vacuum effects 

and matter effects 

are dominating in the past 



-1 - 0.5 0 0.5 1 1.5 2 

X 


Figure 7. The phase portrait for dynamical system (3.11)-(3.13) for real a and H > 0. The grey 
domain represents non-physical solutions. The phase portrait is organized by three critical points: the 
de Sitter universe represented by a stable node (point 1), the Einstein-de Sitter universe represented 
by an unstable node (point 2) and the generalization of Einstein-de Sitter represented by a saddle 
(point 3). 


System (3.11)-(3.13) is three-dimensional but it has invariant submanifolds Y = 0 and 
Z = 0. The behavior of trajectories on the invariant submanifold Y = 0 describes fully the 
global dynamic. The phase portraits on the plane (. X, Z ) are presented in figures 7 and 
8 . Because of the constraint Y 2 = 1 — X — a 2 Z 2 the physical trajectories lie in the region 
Y 2 > 0. Beyond this region is situated a non-physical region (the shaded region in figures 7 
and 8). The boundary of the physical region is determined by a parabola X = 1 — a 2 Z 2 and 
a line X = 0. 
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-1 - 0.5 0 0.5 1 1.5 2 

X 


Figure 8. The phase portrait for dynamical system (3.11)-(3.13) for imaginary a and H > 0. The 
grey domain represents non-physical solutions. The phase portrait is organized by two critical points: 
the de Sitter universe represented by a stable node (point 1) and the Einstein-de Sitter universe 
represented by an unstable node (point 2). 


<*(t) 



Figure 9. Diagram of S(t). The top function represents the evolution of S assuming the best fit of 
H 0 = 68.27 km/(s Mpc) and right boundary values of 95% C.L. for fl mi0 = 0.2416, = 0.0039, the 

middle function represents the best fit for all parameters. The top function represents the evolution 
of 5 assuming the best of H 0 = 68.27 km/(s Mpc) and the left boundary values of 95% C.L. for 
fl m ,o = 0.3542, fi a q = —0.0056. On the t-axis we use a unit times 100 Mpc s/km. See also table 2. 

All critical points lie on this boundary. From the physical point of view they represent 
asymptotic states of system (3.11)-(3.13) which are started at r —> —oo and reach the critical 
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points at t = +oo. The critical point marked as (1) represents the de Sitter universe and 
is a stable node. It is a global attractor for trajectories from its neighborhood. The critical 
point (2) is an unstable node and it represents the CDM model. 

The novelty on the phase portrait is the presence of critical point (3). It is of saddle 
type. A this critical point p m = 2 ^ 1+ ^~— H 2 and Ht = ^ 1 + 3 q2 f he., it represents a 
universe dominated by both decaying vacuum and matter. 

The de Sitter state (critical point 1) is connected by an outcoming separatrix with the 
saddle (point 3). The second separatrix gets in the saddle (point 3) and gets out from the 
Einstein-de Sitter state (point 2). The other trajectories in a non-shaded region start from 
the Einstein-de Sitter state and finish in the de Sitter state. In all cases the time flows from 
r = —oo (a = 0) to r = +oo (a = +oo). At the critical point (3) the decaying A and matter 
play important role and cannot be neglected. At this critical point scale factor a and p m 

2 a 

behaves like: a oc jj = t -i an d p = t - 2 _ 

Vl+oQ;^ — 1 Vl+oQ! z — 1 

This critical point exist only if a / 0. If a — > 0 then it coincides with the CDM 

universe. This critical point represents a new generalized CDM model in which p m = |t -2 
2 

and a[t) oc t 3 in the early universe. 

If the function 


5 = - 


dA(t) 

dt 


2 a 2 
fc 3 


a 


H Pm H p n 


i‘ 3 »o 3 Sf + 


3 y/ n A,0 H 0 

2 


3 y/ n A,0 H 0 


^A ,0 + 


t 2 


1 - 2 n 
3Hot 


+ y/Q 


A.O-7 7== -\— 


2 \ - 1 


3^/f2A,oi?o 

/ J 


(3.15) 


is slowly changing then 

Pm = Pm,oa _3+,5(<) - (3.16) 

Let 5(t) = S = const then p m = p m ,oo _3+<5 - If t —> oo then 5(t) —> const. At the critical point 
(3) 5(t) = 6 = (\/l + 3a 2 — l) 2 /a 2 . The diagram of S(t) is presented in figure 9. 

Some interesting interpretation of our postulated A (t) relation can be derive if we apply 
Starobinsky’s argument [22] that p$ after some averaging over time in the interval At m~ l 
assumes the following form in the quintessence epoch 


p ( / ) = V 0 + Aa 3 


(3.17) 


Therefore in the matter dominating phase we obtain the A (t) parametrization (1.6). 
Finally the model involved belongs to the class of models with so called early dark energy 
constant in which = const = D e during the matter dominated stage (the same refers to 
the radiation dominated stage, too, but with a different value of D e ). 

If 5 <C 1 for the fractional density of dark energy D e [23, 24] it assumes in the interme¬ 
diate domain of the universe the following form 


1 - Q e (a(t)) 
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In the early universe this value is constant 


a. = 


^a,0 

^m,0 ^ 


a,0 


(3.19) 


Therefore for a small value of f2 Q) o, H e is obtained as H e = ^°’° 0 . 

Hojjati et al. [25] found the fraction in total density contributed by early dark energy 
which is approximately equavalent to He¬ 
in our model the exact form of H e (f) is 


H e (f) 


Ha,o H- 
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3Hot 


+ 


In — 1 


A,0" 
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(3.20) 


The evolution of fractional density of dark energy H e has a shape of the “logistic” curve. For 
convenience, the diagram 1 — H e (log(o)) is presented in figure 10. 


l-n e (log(a)) 



Figure 10. Diagram of 1 — H e (log(a)) which is a share of energy density of matter in the total energy 
density. The function represents the evolution of 1 — for H 0 = 68.27 km/(s Mpc), H m> o = 0.35 
and n a>0 = 0.05. The present epoch is at log(a) = 0. For negative values of log(a) we have past 
evolution of 1 — H e (log(a)) with a constant phase of the fractional density of dark energy at the 
early universe (f2 e is small and close to zero). For positive values of log (a) we have future evolution of 
1 — f2 e (log(a)) with a constant phase of the fractional density of dark energy fl e at the late universe 
(H e is big and close to one). Between these two constant phases there is an intermediate phase of 
changing 1 — H e (log(a)) in which we are living. 

There is another interesting approach to running cosmologies, proposed by Starobinsky 
[26] . In this approach (the bottom up), the universe in the quintessence epoch is described by 
a scalar field minimally coupled to gravity with some self-interacting potential. He proposed 
the reconstruction of this potential from the evolution of the scalar perturbation (or also 
luminosity function) in dust like matter component. 

Masso et al. [27] discussed some aspects of contribution to the dark energy density of 
coherent scalar field oscillation in the potential. They obtained using the analytical method 
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of adiabatic invariance that for a quadratic potential the energy density p^ evolves as a^ 3 
and a quartic potential V(</>) ~ evolves like for the radiation matter a -4 . Therefore if we 
add Abare to the potential in a matter (or radiation) dominating universe p$ = A(i). 


4 Statistical analysis of the model 


In this section we present a statistical analysis of the model parameters using the SNIa, BAO, 
CMB observations, measurements of H(z) and the Alcock-Paczynski test. 

First, we use the Union 2.1 sample of 580 supernovae [28]. For the SNIa data we have 
the following likelihood function 


In UsNIa 



(4.1) 


where the summing is over the SNIa sample; the distance modulus jU obs = m — M (where 
m is the apparent magnitude and M is the absolute magnitude of SNIa stars) and p th = 
5 log 10 Dl + 25 (where the luminosity distance is Dp = c(l + z) J 0 " jj^ and a is the uncer¬ 
tainties. 

We use the BAO (baryon acoustic oscillation) data which were taken from the Sloan 
Digital Sky Survey Release 7 (SDSS DR7) dataset which consists of 893 319 galaxies [29]. 
The likelihood function is given by 


In Tbao 
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r s (~d) 

D v (z) 


d{z) j 
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(4.2) 


where r s (zd ) is the sound horizon at the drag epoch and z = 0.275, d(z) = 0.1390, a = 0.0037 

[ 3 °]- 

The next likelihood function encompasses the Planck observations of cosmic microwave 
background (CMB) radiation [31], the information on lensing from the Planck and low-£ 
polarization from the WMAP and has the form 


hlLcMB+lensing+WP = ~ V* - X° bs ), (4.3) 

ij 


where C is the covariance matrix with the errors, x is a vector of the acoustic scale l a, the 
shift parameter R and f Ibh? where 


7 r f dzf 

A ~^ r ) c J 0 ~W~Y 


(4.4) 

(4.5) 


where z* is the recombination redshift and r s is the sound horizon. 

The idea of the Alcock-Paczynski test is the comparison of the radial and tangential 
size of an object, which is isotropic in the correct choice of model [32, 33]. The likelihood 
function is independent of the parameter Hq and has the following form 


In L A p 



(AP th ( Zi ) - AP ohs (zi))' 2 


(4.6) 
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where AP(z) th = / 0 " and AP(zi) ohs are observational data [34-42]. 

At the end it is also valuable to add the constraints on the Hubble parameter, i.e. 
H(z = 0) = H 0 . 

Data of H{z) for samples of different galaxies were also used [43-45]. 


In L H ( z ) 



(4.7) 


The final likelihood function for the observational Hubble function is 


Ltot = TsNIa-hBAO-^CMB+lensing+WP^ApAif( z ). 


(4.8) 


^tot 



Figure 11. Diagram of PDF for parameter Hq in units km/(100 s Mpc) obtained as an intersection 
of a likelihood function. Two planes of intersection likelihood function are = 0.2938 and n a ,o = 
-0.0006. 

To estimate the model parameters we use our own code CosmoDarkBox implementing 
the Metropolis-Hastings algorithm [46], [47]. We use the dynamical system formulation of 
model to obtain the likelihood function [30, 48]. 

We use observation data of 580 supernovae of type la, selected subsets of the data points 
of Hubble function, the measurements of BAO from SDSS DR7. We also use data for the 
application of the Alcock-Paczynski test 18 observational points. At last, we estimated model 
parameters with CMB data from Planck, low- l polarization from WMAP and lensing from 
Planck. To estimate the model parameters we chose interval (64.00, 74.00) for Hq and (0.21, 
0.37) for D m! o- The values of estimated parameters for a 2 from the interval (-0.05, 0.05) are 
shown in table 2 and for positive a 2 from the interval (0.00, 0.05) are shown in table 3. The 
best fit for model with a 2 from the interval (-0.05, 0.05) is in the part of likelihood function 
where a 2 is negative. 

If it is chosen the lower limit of the interval of a 2 larger than the value of the best fit 
then the best fit of the model for the new interval is equal the value of the lower limit of this 
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^tot 



Figure 12. Diagram of PDF for parameter o obtained as an intersection of a likelihood function. 
Two planes of intersection likelihood function are fi mj o = 0.2938 and Hq = 68.27 km/(s Mpc). 


^tot 



Figure 13. Diagram of PDF for parameter fl m! o obtained as an intersection of a likelihood function. 
Two planes of intersection likelihood function are fl aj o = —0.0006 and Hq = 68.27 km/(s Mpc). 


interval. A consequence is the lower limit of the error is equal zero for a 2 . So the specific 
values of the best fit and errors of f2 ai 0) a 2 and jo in table 3 are a result of the choice of 
limits of the interval of a 2 . 

To illustrate the results of statistical analysis the diagrams of PDF are shown in fig¬ 
ures 11, 12 and 13. In turn figures 14 and 15 shown the likelihood function with 68% and 95% 
confidence level projection on the (D a; 0 ) ^m,o) plane and the (D 0i o, Ho) plane, respectively. 
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Figure 14. The likelihood function of two model parameters (f2 a> o, fl m ,o) with the marked 68% and 
95% confidence levels. The value of Hubble constant is estimated from the data as best fit value 
Hq = 68.27 km/(s Mpc) and then the diagram of likelihood function is obtained for this value. 


Table 2. The best fit and errors for the estimated model with a 2 from the interval (—0.05, 0.05). 


parameter 

best fit 

68% CL 

95% CL 

H 0 

68.27 km/(s Mpc) 

+0.96 

-1.07 

+1.67 

-1.68 

0 

0.2938 

+0.0355 

-0.0325 

+0.0604 

-0.0522 

^O!,0 

-0.0006 

+0.0031 

-0.0030 

+0.0045 

-0.0050 

jo 

1.002 

+0.010 

+0.016 

-0.008 

-0.013 

a 2 

-0.002 

+0.010 

-0.007 

+0.014 

-0.012 


5 Conclusion 

The aim of the paper was to study the dynamics of the emerging A(f)CDM cosmological 
models. In the study of dynamics we find exact solutions and use dynamical system methods 
for the analysis of dependence of solutions on initial conditions. In the latter evolutional paths 
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Figure 15. The likelihood function of two model parameters with the marked 68% and 

95% confidence levels. The value of fl m ,o constant is estimated from the data as best fit value 
fl m ,o = 0.2938 and then the diagram of likelihood function is obtained for this value. 

Table 3. The best fit and errors for the estimated model with positive a 2 from the interval (0.00,0.05). 


parameter 

best fit 

68% CL 

95% CL 

Ho 

68.38 km/(s Mpc) 

+0.97 

-0.97 

+1.16 

-1.16 

^m,0 

0.2877 

+0.0198 

-0.0242 

+0.0331 

-0.0442 

^O!,0 

0.0000 

+0.0025 

-0.0000 

+0.0038 

-0.0000 

jo 

1.000 

+0.000 

+0.000 

-0.008 

-0.011 

a 2 

0.000 

+0.008 

-0.000 

+0.012 

-0.000 


of cosmological model are represented by trajectories in the phase space. Due to geometrical 
visualization of dynamics we have the space of all solutions and can discuss their stability. 
We are looking for such trajectories for which the ACDM model is a global attractor in the 
phase space. 


19 - 



















We study in details dynamics of cosmological model with the prescribed form of A(t) = 
Abare+yhu where Abare is a positive constant and a 2 is either positive or negative. We calculate 
exact solutions for the scale factor and subsequently calculate the jerk. It is demonstrated 
that this parameter is time dependent if and only if the effects of time contribution to A (t) 
are non-zero. We propose the measurement of the jerk as a diagnostic of decaying A, i.e. 
A < 0. Due to analysis of dynamics in the phase space we have found an interesting solution 
in the phase space, a saddle critical point, at which p m (t) scales like t~ 2 . This solution was 
recently proposed by Ringermacher and Mead [16] as a description characteristic for the dark 
matter evolution. 

From the phase portrait we derive the generic scenario for an evolution of cosmological 
models with such a form of the dark energy parametrization. Trajectories start from initial 
singularity (the Einstein-de Sitter model) and then go in vicinity of the saddle point where 
they spend a lot of time and then go to the de Sitter state. It is a typical behavior for all 
generic trajectories in the phase space. The new critical point is emerging in the phase space 
due to the effect of a time dependence of the cosmological constant. If a 2 < 0 then this point 
is absent (it is gluing with the critical point representing the Einstein-de Sitter model). 

We also tested this model using astronomical data. Statistical estimations show that 
the model fits to data as well as the standard cosmological model (the ACDM model). In any 
case the value of a 2 = 0 belongs to the confidence interval for the estimated parameter a 2 we 
cannot reject that a 2 ^ 0. Only if we find the best fit value a 2 with the error of the one order 
less than this value the problem of estimated / cxq / 0 could be solvable. We can obtain the 
limits on the value of parameter a 2 and —0.009 < a 2 = 3HqTq (1 — D m) o — Da,o) < 0.008 (for 
68 % C.L.) and —0.014 < a 2 < 0.012 (for 95% C.L.). 

In papers of Doran and Robbers [23] and Pettorino et al. [24] there are limits on 
fractional dark energy at early time. Recently Ade et al. [49] have found D e < 0.0036. It is 
interesting that they have obtained a similar limit to our limit on Q, a < 0.0038 in the other 
parametrization of dark energy. 

Note that if we apply Starobinsky’s idea the parameter a 2 can be constrained through 
D e measurement. This parameter measures amount of dark energy at the early evolution 
of the Universe. If D e is different from zero then we obtain value information about this 
alternative evolutional scenarios which are consistent with the present epoch. 

In our case D e = < 0.0036 and therefore a 2 < 3 HqTq D m .oD e . If we put 

D m ,o = 0.25 and HqTq = 1 then we obtain a 2 < |O e = 0.0027. 

Finally we obtain a stronger limit for a 2 then in table 3. However, note that this esti¬ 
mation is model dependent (it is assumed Starobinsky’s argument), note that the case study 
of our model fully confirm existence of phase during the early universe at which fractional 
energy density of dark energy is constant (see figure 10 and eq. 3.20). 
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